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Suppose (hat 

P n = 10 log 10 [10 Vl/, ° + ••• + 10- W10 ], 

where \X n \ is a sequence oj independent random variables. The main 
result of this paper shows that under very general conditions on the sequence 
{X n \, the power sums P n will be asymptotically normally distributed. 
This result supports a commonly used normal approximation, and shows 
why many physical quantities obtained by power addition oj random variables 
tend to be normally distributed in dB. 

I. INTRODUCTION 

In many areas of transmission engineering, logarithms of sums of 
powers are considered in the form 

P n = 10 log 10 [10 Vl/10 + ••• + 10 W10 ], 

where X lt . . ., X n are random variables. Specifically, if X\, . . ., X n 
are power levels in dB such that 

X, = 10 log.o (Wj/w,,) j = 1,2, ■■■ ,n, 

where w Q , w x , ■ ■ • , w„ are powers (e.g., expressed in watts), then the 
power level in dB of the sum w = w x + • ■ • + w n is given by the so- 
called "power sum," 

P n = 10 log 10 (w/uO = 10 log 10 [10 Y ' /10 + • • • + 10 A '" /10 ]. 

Quite often X u • • • , X„ are taken to be mutually independent random 
variables with specified distributions, and it is of interest to determine 
properties of their power sum P„. 

A major difficulty encountered in working with power sums is that 
the distribution and moments of such a sum usually cannot be ex- 
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pressed in simple closed form. This includes, for example, the im- 
portant case when X i} ■ ■ • , X n are mutually independent and each has 
a truncated normal distribution. Even in the simpler case when X 1} 
• • • , X„ are mutually independent, identically distributed, and Xi is 
normal, the problem is intractable. The difficulty and importance of 
the general problem, in turn, has led to a number of methods for ap- 
proximating the distribution of a power sum. 1 - 2 - 3 - 4 - "• °- T - 8 - 9 

In the present paper, the asymptotic distribution of a power sum 
is studied. The main result is a limit theorem which shows that under 
very general conditions on the components X 1} X 2 , • • • , the correspond- 
ing power sums P„ will be asymptotically normal as n -> oo. The par- 
ticular form of the result is as follows: Given a sequence {X„) of 
mutually independent random variables satisfying certain conditions, 
there exist sequences of constants {c„} and {(/„} such that 

limP{[(P„ - c n )/d n ] ^ x] = [1/V^] [ exp [-f/2] dt. (1) 

The conditions for (1) to hold are the central concern of this paper, 
but the implications of the results are equally important. In particular, 
one of the oldest and most useful approximations to the distribution 
of a power sum is a normal approximation. This approximation was 
first used at Bell Telephone Laboratories in 1934 by R. I. Wilkinson, 2 
and is based on the fact that many observed power sum distributions 
are "nearly normal." This includes power sum distributions obtained 
by numerical convolution, and empirical distributions of physical 
quantities such as noise levels on trunks and connections where the 
resultant noise (on a dB scale) can be viewed as an approximate power 
sum. 10, J1 The limit theorem proved in this paper thus provides mathe- 
matical support for a normal approximation, and substantially explains 
why many physical quantities obtained by power addition of random 
variables tend to be normally distributed in dB. 

II. A NORMAL LIMIT THEOREM FOR POWER SUMS 

2.1 Discussion 

Before stating the main results, it is instructive to show informally 
why one would expect power sums to be asymptotically normal. To 
take a simple case, suppose that {X„} is a sequence of mutually inde- 
pendent, identically distributed random variables such that 

r 2 - Var [10- Yl/1 °] 
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is finite. Let 6 = E10 x,/l ° and put 

S n = 10 v,/l " + • • • + 10 Vu/1 °. 
Then by the law of large numbers, one expects that for large n, 

Next, note that if .<• ;^ 1, then log. x Rrf x — 1 so for large n 

, S„ S n - nd 

lo g' : n ~ 



// e n e 

Multiplication by (0 Vn ) ..V then gives 



eVn , S n _ S n - nO ... 

t »0 rVn 



But, by the central limit theorem, the right-hand side of (2) is asymp- 
totically normal with mean and variance 1 . Thus, it is strongly sug- 
gested that 

lim p\^y~^ [log, S n - log, (nd)] ^ x 

= [l/Var] f exp[-t 2 /2]dt. 
This, and more, is indeed true as will be shown. 

2.2 The Main Result 

The normal limit theorem for power sums is a consequence of the 
following result which will first be proved: 

Lemma 1: Let {S n } be a sequence of positive random variables. Sup- 
pose there exist sequences of positive real numbers {a,,} and {b,,}, and 
a distribution F such that 

(i) At each point of continuity of F, 

(«") lim (b n /a n ) = 0. 

Then at each point of continuity of F, 

limPJ(fl n /6 n ) log, (S n /a n ) £ x) = F(x). 
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Proof: Let z be a continuity point of F, and let e > be given. Because 
F has at most a countable number of discontinuities, there is a 8 > 
such that F is continuous at x + 8 and 

F(x + 8) - F(x) < c. (3) 

Next, define 

U n = (& - a n )/b n and V n = (a n /b n ) log, (SJa H ). 
Then 

\P[V n SZ) - F(X) | 

^ |P{7„ ^ x] -P{U n £x) |+ |P{tf» ^ z) - F(s) |. 
By assumption (i) therefore, 

lim"|P{7 n ^z} - F(x)\ ^ Urn |P{7„ ^ x) - P{tf„ ^ *} |. 

Let 

A„(x) = |P{F n ^z} - P[U n ^ x) |. 
To complete the proof it suffices to show that 

lim A„(x) = 0. 

n-»oo 

To prove this note first from the inequality log. x g X - I, x > 0, that 
V„ g U n for all n. Thus, 

a„(.t) =p[{f b ^ini^> »j] 

= P{x < U n ^ (a n /6 n )[exp (b n x/a„) - l)\ . 

Using the inequality e v - 1 ^ 2/e", - « < 2/ < » , it follows that 

^ A n (.-c) ^ P{x < U n ^ x exp (b n x/a n )}. 

By assumption, (6„/a n ) > for all n and lim,,-*, (b n /a„) = 0. Thus, 
there exists a natural number N such that n ^ N implies 

x < x exp (b n x/a n ) ^ z + 5. 

So if n ^ iV, 

^ A n (.r) ^ P{x < U n ^ x- + 5j . 

Because a; and x + 8 are continuity points of F, it follows by assump- 
tion (i) and inequality (3) that 
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g Tim A„(.r) ^ F{x + 5) - F(x) < e. 

n— co 

Since e > was arbitrary, the proof is complete. 

The importance of Lemma 1 is that it gives a sufficient condition to 
go from limit theorems for sums of random variables to limit theorems 
for logarithms of sums. In the important case of power sums of in- 
dependent random variables, general conditions for asymptotic nor- 
mality can thus be obtained from classical central limit theory as 
shown in the next result. 

Theorem 1: Let {X n } be a sequence of mutually independent random 
variables and suppose that 

t) s Var [10 Y,/10 ] 
is finite jor every j. Let 0, = E10 X ' /10 and put . 

M n =te,, si =±t% 

i-l j'-l 

Denote the distribution of 10 A,/1 ° by Hj(x), and let 

P n = 10 log 10 [10 v,/, ° + • • • + 10 A '" /10 ]. 
// the following conditions are satisfied: 

(i) The Lindeberg Condition: For every e > 0, 

lim I, £ f (x - e,) 2 dH f (x) = 0, 

where 

A in = [x: | x — 6j | ^ es n ) 

(ii) lim (s n /M n ) = 

it will follow that 

lim P{(XM n /s„)[P„ - 10 log 10 M n ] ^ x] = *(.r) (4) 

n—co 

where A = (log c 10)/10 and 

*(.r) - [1/V^] f exp [-f/2] dt. 

J — 00 

Proof: Let 

S n - 10 Xl/1 "+ ••• + io- W10 . 
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Then condition (i) implies that 

.. JS n - M n , \ „,, 

hmP< — ^ x> = *(.t) 

(cf. Feller, 12 p. 256). With the identifications a„ = M n and b n = s„ it 
follows from condition (u) and Lemma 1 that 

limP{(M„/s n ) log. (SJM n ) ^ x} = $(*). 

n-»oo 

The assertion of the theorem then follows by changing to logarithms 

with base 10. 

An interesting thing to note is that if the conditions of Theorem 1 
are satisfied then the sum of powers 

S n = 10 x,/1 ° + • • • + 10*" /10 
and the power sum in dB, P„ = 10 log 10 <S„, will both be asymptotically 
normal. Thus, not only will normality be observed on a "power scale" 
but on a "dB scale" as well. 

2.3 Identically Distributed Components 

The preceding result implies the asymptotic normality of P„ when 
the components are identically distributed. To show this, suppose that 
{X n \ is a sequence of mutually independent, identically distributed 
random variables with H (x) = P{10* ,/10 g re}. Let 

r 2 = Var [10- v,/1 °] 

and 6 = E10 x ' no . If t 2 is finite, condition (u") of Theorem 1 is clearly 
satisfied since 



M n eVn 
Condition (i) is also satisfied because if e > 0, 

ht, f (x- e,) 9 dH,(x) = \ [ (x - &? dH(x) -> as n -> «, , 

S„ fTl J Aln T J An 

where A n = [x : | x — | ^ er y/n }. It thus follows that 

lim p{\ ^ [P, - 10 log 10 (nd)] ^ *} = -^= £ exp [-f/2] dt, 

hence, P„ is asymptotically normal with mean 10 log io (n0) and vari- 
ance T 2 /{n\'-6-). 
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2.4 Bounded Components 

Suppose next that {X,,} is a sequence of mutually independent ran- 
dom variables and that the following conditions are satisfied: 

(i) There exist constants b and B such that 

< b ^ 10 v,/1 ° ^ B for all j 

(ii) si — » oo as n — » °° . 

The conditions of Theorem 1 arc easily shown to be satisfied in this 
case, and it follows that P n will be asymptotically normal. Note that 
condition (i) will be satisfied whenever 10 A,/1 ° represents power from 
a physical source. Condition (/•/), on the other hand, will be satisfied 
if t 2 ^ c > for some fixed c and an infinite number of indices j. 

III. THE NORMAL LIMIT THEOREM AND WILKINSON'S 
NORMAL APPROXIMATION 

One of the most useful approximations to the distribution and mo- 
ments of a power sum is based on a normal approximation as men- 
tioned in the introduction. The method consists of approximating the 
distribution of P„ by a normal distribution so that 

P{1\ < x) «P[^ T & ^ x}, 
where $ is normal with mean and variance 1. Writing as before, 

M n = y^lO 7 '" 710 and si = Var [10 p " /10 ], 
the parameters a and ft are chosen so that 

M n = £[10 (a ^ ,/10 ] 
and 

s 2 n = Var[10 ( " £+3,/1 °] 

which is equivalent to equating means and variances on a "power 
scale." If £ is normal with mean and variance 1 then 

JS?[10 (a * + " )/l0 ] = e x V (Xa) ' 
and 

Var [10 ( "« +/,)/1 °] = c^[e 2X "" -e x ""], 
where 

X = <>g. 10)/10. 
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Solving the above equations for a and (3, the approximation then as- 
serts that P n is normal with 

E(P n ) = p = 10 log 10 M n - 5 log 10 [1 + (s n /M n ) 2 ] (5) 

and 

Var CPJ - a = f log 10 [1 + (s R /M n ) 2 ]. (0) 

In light of the normal limit theorem, it is quite natural to assume 
that P„ is approximately normal, provided the conditions of the 
theorem are satisfied, and n is large. On the other hand, the estimates 
given by (5) and (6) are different from those based on (4) : 

E(P n ) = 10 log 10 M n (7) 

Var(P n ) =sl/(\M n )\ (8) 

The difference, however, is easily resolved once it is realized that if 
condition (tt) of Theorem 1 is satisfied then (5) and (6) are asymptot- 
ically equivalent to (7) and (8). In fact, it is a simple matter to 
show (cf. Feller, 12 p. 246) that if the conditions of Theorem 1 arc 
satisfied then 

limP{[(2\ - u„)/V7 n ] ^ x\ = -if* f exp [-i 2 /2] dt, (9) 

n-oo v Zir J -<=o 

where 

u n = 10 log 10 M n - 5 log 10 [1 + (sJMS'] 



and 



v. = f log 10 [1 + (s n /M n )*]. 



In numerical applications, the normal approximation based on (9) 
is to be favored over that based on (4). In the first place, when 
Xi , • • • , X„ are mutually independent, identically distributed, and A r x 
has a truncated normal distribution, Monte Carlo studies by I. Nasell 
have shown that the mean and variance estimates given by (5) and (6) 
are better than those given by (7) and (8) (although for large n and 
small variance of A'i there is hardly any difference). Secondly, the nor- 
malizing factors in (9) were obtained quite naturally by equating mo- 
ments on a power scale. This is analogous to the situation in classical 
central limit theory when the sequence (S n —M„)/s„ converges in dis- 
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tribution to the standard normal. The normalizing factors M„ and s„ 
are not the only ones that give this result, but they are chosen in a 
natural way to insure that for every n, the mean and variance of 
(S n —M n )/s n agrees with its asymptotic distribution. 
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